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Abstract 

We study spin-dependent transport through a spin diode in the pres- 
ence of spin-flip by means of reduced density matrix approach. The cur- 
rent polarization and the spin accumulation are computed and influence of 
spin-flip on the current polarization is also analyzed. Analytical relations 
for the current polarization and the spin accumulation are obtained as a 
function of polarization of ferromagnetic lead and the spin-flip rate. It is 
observed that the current polarization becomes zero under fast spin-flip 
and the spin accumulation decreases up to %85 when the time of spin-flip 
is equal to the tunneling time. It is also observed that the current po- 
larization increases linearly when the dot is singly occupied, whereas its 
behavior is more complicated when the dot is doubly occupied. 

1 Introduction 

Spin dependent transport through mesoscopic systems has been widely studied 
both experimentally and theoretically points of view [TJ HJ [3J HI [5] . Study of 
transport through quantum dots (QDs) provides worthful information about 
novel physical phenomena such as Kondo effect [BJ [7] , spin and Coulomb block- 
ade effects [5J |H1 [13 HI] ) tunneling magnetoresistance [HJ[T2J[13], etc. Coupling 
of the QD to different leads can result in different transport properties. Spin 
filters composed of normal metallic leads [TS1 [H] and ferromagnetic leads have 
been applied for producing spin accumulation affecting on negative differential 
resistance (NDR) [T7] and zero bias anomaly |18) . 

An interesting structure is composed of the QD coupled to a nonmagnetic 
lead and a ferromagnetic lead. It has been shown that the QD can act as both a 
spin diode and a spin filter by reversing bias voltage |14[|17j . Diode- like behavior 
in a carbon nanotube coupled to a normal metal and a ferromagnetic lead was 
recently reported [19.. This behavior was also observed in resonance tunneling 
connected to ferromagnetic contacts with different spin-dependent transparen- 
cies |20j . Spin diodes have been theoretically investigated without taking into 
account spin-flip [211 22;. In this paper, we analyze this structure with taking 
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into account spin-flip process in sequential tunneling regime. Rate equations are 
used to describe the states inside the QD with considering coherency between 
different states. This subject was recently studied by C. Feng and co-workers 
[2"5] . Here, we extract analytical relations for spin current polarization and spin 
accumulation. We show how spin-flip decreases the current polarization and 
its effect on the spin accumulation is also studied. Dependence of the current 
polarization on the polarization of the magnetic lead is also examined. 

The article is organized as follows: in the next section, we compute the 
spin-dependent current using rate equations. Then, the numerical results are 
presented and the behavior of the system is analytically studied. Finally, some 
sentences are given as a conclusion. 

2 Description of the model 

We consider a single level quantum dot coupled to a normal metal and a ferro- 
magnetic lead. Hamiltonian describing the system is 

H = Hi ea d + Hd + H s .f + Ht (1) 

where Hi eads = J2 a ka £ akaC i akcr c ak(T describes the leads, c aka {c\ ka ) destroys 
(creates) an electron with wave vector k, spin a in the lead a (a — L,K). In 
continue, it is assumed that the left lead is normal metal while the right one 
is ferromagnetic. Hence energy of electron is spin-independent in the left lead, 
whereas for the right lead e aka = e ka — (l) s "->-A that A is band spin splitting. 
Hd = Ylcr £ d n <j + Un^n^ where U is Coulomb repulsion, and n a = d\d a is 
occupation operator. d (7 (dJ T ) denotes the annihilation (creation) operator in the 
QD. Ed = £o — cugVg — ollVl — aaVn is the energy level of the QD. We assume 
that the QD is capacitively coupled to the gate (G), left and right electrodes. 
Ht = J2ak<TiVaLkaC aka d a + H.C] describes the tunneling between leads and the 
QD. H s j contains spin-flip scattering inside the QD, H s j = R[d].di + d/^df] 
where R is spin-flip rate. This process can be obtained by, for instance, spin- 
orbit interaction [53] or a transverse magnetic field [T5] that rotates the electron 
spin. 

The QD can be in i) empty state |0 >, ii) singly occupied state \a > (a =t, I) 
( £ \a> = £ d) or iii) doubly occupied state \2 > (e|2> = 2ed + U). It is clear that 
the states | t> an d | i> are not the eigenstates of the isolated QD because 
of spin-flip. In order to study the QD, reduced density matrix has been used. 
Using Markov approximation, time evolution of reduced density matrix elements 
are given as follows [551 (fi = 1) : 

-i < s\[H s . f ,P]\s' > +6 SS , J2( W ksPkk ~ W sk P ss ) (2) 

k^s 

i[i - Wsk + w s>k]p ss > 



dP ss > 
dt 
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P ss is the probability of being in the state s (s = 0, a, 2), whereas P ss > describes 
the coherency between the states \s > and \s' >. [A,B] = AB — BA and Wks 
denotes the transition from \k > to \s > and is computed by means of Fermi's 
golden rule as 

W ks = ^2 T -lU\ s ks\)S Nk ,N s+ i + f- (|e*.|)<W.-i] (3) 

a 

where stands for the number of electrons in the state \k > (Nk — 0, 1,2) 
and Eks — £k — £s is the transition energy. f a (x) = [1 + exp((x — ji a ) / 'fcT)] _1 is 
Fermi distribution function of the a th lead with chemical potential a, whereas 
/~ = l—f a . For computing the tunneling rates (Wks), wide band approximation 
is used i,.e. dot-lead coupling is energy independent and = ^Z k \V a ka\ 2 ■ 
Hence, for instance, Poo and Pfi are obtained from eq.2 as 

Y^W^Paa-WoM (4a) 
a 

-iR[P u - P n ] -\[J2 W ^ + J2 W ik\Pn (4b) 

fe/t kjil 

and it is clear that P^ = P^. 

Solving Eqs. 2 in the steady state (dP/dt — 0), spin-dependent current 
crossing from the lead a is computed by 

la - r-[/ Q ( £d )P 00 - f~P aa + f a {e 2s )P ss - f-(e 2s )P 22 ] (5) 

where a is the opposite of a. For simulation purpose we set T a L = To and 
F^ = (1 + (— l) s ^- a p)T where r is the dot- lead coupling strength and p is the 
spin polarization degree of the right lead. We also set /il =0 and fiR = — V 
where V is applied bias. Therefore, the left lead acts as an emitter when the 
bias is positive, while the right one acts as the emitter in negative bias. 



dPpo 
dt 

dPu 
dt 



3 Results and discussions 

Fig. la shows the current polarization £ = jrifji as a function of bias. As 
we expect, spin-flip reduces the current polarization because this process tries 
to destroy the spin polarity of the system. It is observed that the device can 
operate as a spin diode when the QD is singly occupied. In positive bias when 
< £d < f-L, the current polarization becomes zero because the left lead acts 
as the emitter. Note that I a = T^poo and for the left lead, we have r£ = rj^ 
and as a result J 1 * = P. On the other hand, for negative bias i.e. \il < £d < Vr, 
the current polarization is maximum because of T^ R > rjj. Indeed, a spin- up 
electron enters the QD faster than a spin-down electron from the right lead. 
Although the device acts as a spin diode in the presence of spin-flip, the current 
polarization reduces significantly. When the QD is singly occupied, the current 
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polarization is given as [27J 



Aa 2 



1 



/or V > 
/or V < 



(6a) 
(6b) 



where a = R/Tq- It is clear from above equation (Eq. (|6bjl ) that increasing a 
gives rise to decreasing £. Note, for a = 0, our results are the same as Eqs. 21 in 
the Ref. [3T] . If the time of spin- flip is much shorter than the time of tunneling 
i.e. a —> oo, the current polarization becomes zero. Indeed, spin-flip destroys 
the effect of the ferromagnetic lead on the QD, completely. 
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Figure 1: (a) The current polarization versus the bias voltage for R = (solid), 
R = T (dashed) and R = 2T (dash-dotted). Parameters are Vg = lmV, 
U = AmV, a G = 1, a L = a R = 0.5, T = 2.5K, T Q = 20^eV and p = 0.5. The 
current polarization as a function of the right lead polarization when (b) the 
QD is singly occupied and (c) the QD is doubly occupied. V = 6mV (solid) 
and V — —6mV (dashed). We set R — (blue) and R = T (gray). 



The dependence of the current polarization on the p under conditions that 
the QD is singly occupied is plotted in the fig. (lb). In positive bias, the 
current polarization is always zero, except in p = 1. When the right lead is a 
half metal the current polarization will be equal to one because the spin-down 
electron injected from the left lead can not leave the QD, therefore 1^ = 0. This 
behavior was reported before in the Ref. [21] . It is straightforward to show that 
the probability of being in the state \a > is given as [27] : 

4a 2 + 1 + (-lj^tp 

3[l+4a 2 ]-p 2 {> 
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It is clear from the Eq. under conditions that R = 0, the QD is occupied 
by a spin-down electron. In the presence of spin-flip, this electron can change 
to a spin- up electron. More specifically, if spin- flip occurs fast enough (a > 1), 
we have Pff = Pff. As a consequence, the current polarization is zero under 
fast spin-flip even though the right lead is a half metal. In negative bias, the 
current polarization increases linearly, as we expect from Eq. 5b. Increasing 
p means that the spin- up electron enters the QD faster and , hence If > P. 
The dependence of the current polarization on the p when the QD is doubly 
occupied is plotted in the fig. lc. It is observed that the current polarization 
is the same for both positive and negative biases. It can be shown that in high 
bias regime [28) . the current polarization is 

^ = V^F? (8) 

It is observed that when a = and p = 1, the current polarization approaches 
1. 
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Figure 2: (a) The spin accumulation as a function of bias for R = (solid), 
R = Tq (dashed) and R = 2Tq (dash-dotted), (b) and (c) show the spin 
accumulation when the QD is singly occupied and doubly occupied, respectively. 
V = 6mV (solid) and V = -6mV (dashed). We set R = (blue) and R = T 
(gray). Other parameters are the same as fig. 1. 

The spin accumulation (to = Pff — P^) is shown in the fig. 2a. In positive 
bias, the spin accumulation is negative because the spin-up electron injected 
from the left lead leaves the QD faster and, hence Pff > Pff. In negative bias, 
the right lead operates as the emitter so that the spin-up electron enters the 
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QD faster and as a result the spin accumulation becomes positive. Spin-flip 
scattering reduces the spin accumulation significantly. In the presence of spin- 
flip, it is probable that the electron in the QD rotates. This rotation leads to 
the reduction of the spin accumulation. There are tow plateaus in each side of 
the figure. The first one is due to entering Ed inside the bias window whereas 
the next one is created when 2ed + U takes place inside the bias window. 

The dependence of the spin accumulation on the p is plotted in the figs. 
2a and 2b. It is observed that when the QD is singly occupied and the right 
lead is a half metal, the spin accumulation will be equal to — 1 if no spin-flip is 
considered. This behavior was predictable from Eq. Q. If the spin-flip rate is 
equal to To, the spin accumulation decreases up to %85. Unlike positive bias, 
the spin accumulation increases linearly in terms of p in negative bias. It can 
be shown that in single electron state m is given by 

' for V > (9a) 



p 2 - 3[4a 2 + 1 
2p 



12a 2 



for V < (9b) 



It is interesting to note that in positive bias, even though it is assumed that the 
right lead is a half metal and no spin flip is considered, m is always lesser than 
1. It results from this fact that the probability of finding the QD in the empty 
state is equal to 1/3. Indeed, the probability of being in the state \a > is given 
as: 

_ l+4 Q 2 + (-l)^p 

3[l + 4a 2 ] liUJ 

It is clear from above equation that if we set p — 1 and a — 0, will be 2/3. 
The magnitude of m is the same for both positive and negative bias when the 
QD is doubly occupied. In this situation, one can obtain [28] 

m = T A - A (11) 

4a 2 - p 2 + 4 y ' 

where the minus sign is for positive bias and plus sign is for negative bias. It 
is also found that m is always lesser that 1 because of interplay between spin 
accumulation and Coulomb interaction. 



4 Conclusion 

In this article, we analyze spin current polarization through a QD coupled to a 
normal metal and a ferromagnetic lead by means of rate equations in the pres- 
ence of spin-flip scattering. It is observed when an energy level of the QD is 
inside the bias window, this system can operate as a rectifier for spin current po- 
larization and show how spin-flip can decrease the magnitude of the polarization. 
It is observed that the spin accumulation and the current polarization decrease 
significantly when spin-flip occurs faster than the tunneling event. Dependence 
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of the current polarization and the spin accumulation on the polarization of 
magnetic lead is also studied and analytical relations are obtained as a function 
of the polarization and spin-flip rate. 
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